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Abstract. We explore quantum signatures of classical chaos by studying the rate of information 
gain in quantum tomography. The tomographic record consists of a time series of expectation 
values of a Flermitian operator evolving under application of the Floquet operator of a quantum 
map that possesses (or lacks) time reversal symmetry. We find that the rate of information gain, and 
hence the fidelity of quantum state reconstruction, depends on the symmetry class of the quantum 
map involved. Moreover, we find an increase in information gain and hence higher reconstruction 
fidelities when the Floquet maps employed increase in chaoticity. We make predictions for the 
information gain and show that these results are well described by random matrix theory in the 
fully chaotic regime. We derive analytical expressions for bounds on information gain using random 
matrix theory for different class of maps and show that these bounds are realized by fully chaotic 
quantum systems. 
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1. Introduction: Classical and Quantum Chaos 

Classical chaos is characterized by the sensitive dependence on initial conditions in a 
deterministic dynamical system m. In a conservative Hamiltonian system, this occurs 
for trajectories which do not settle down to fixed points, periodic orbits, or quasi-periodic 
orbits in the limit t ^ oo, where t is the time of evolution of the trajectory lEIIl. Sensitive 
dependence on initial conditions means that nearby trajectories separate exponentially 
fast; the rate of separation is given by the Lyapunov exponent. A, which characterizes the 
dynamics of the system. A conservative Hamiltonian system with N degrees of freedom, 
with N constants of motion, it is said to be integrable, and its dynamics is regular. When 
there are fewer than N constants of motion, then the individual trajectories can explore 
the phase space in a complex manner and the system can exhibit chaos. 

It is not difficult to see that the above “definition” of chaos fails in the quantum do¬ 
main. A quantum state is not a point in the phase space but is described by a state vector. 
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The time evolution of the state vector, due to the Schrodinger’s equation, is unitary. This 
means that the overlap of two state vectors undergoing evolution is constant with time. 
Therefore, quantum systems, unlike their classically chaotic counterparts, do not show a 
sensitive dependence on initial conditions. Furthermore, while classical chaos can lead to 
infinitely fine structures in the phase space, in quantum mechanics, Planck’s constant, h, 
sets the scale for such structures, according to Heisenberg’s uncertainty principle. This is 
often stated as the key reason for the absence of chaos in the quantum domain. This, how¬ 
ever, is not the complete story. An alternate description of classical mechanics, involving 
the evolution of classical probability densities, preserves the distance between two prob¬ 
ability densities as a function of time 0. Hence, the distance between two probability 
densities does not show exponential sensitivity even for classical mechanics. 

All this leads to two interesting questions: 

1. How does classically chaotic dynamics inform us about certain properties of quan¬ 
tum systems, e.g., the energy spectrum, nature of eigenstates, correlation functions, 
and more recently, entanglement and quantum discord. Alternatively, what features 
of quantum systems arise due to the fact that their classical description is chaotic? 

2. Since all systems are fundamentally quantum mechanical, how does classical chaos, 
with trajectories sensitive to initial conditions, arise out of the underlying quantum 
equations of motion? 

These two questions are not unrelated. However, the first question deals mainly with find¬ 
ing the signatures of chaos by studying the properties of the quantum Hamiltonian, while 
the second concerns with the dynamical behaviour of quantum states and the emergence 
of classically chaotic behaviour in the macroscopic limit. 

A central result of quantum chaos is its relationship to the theory of random matrices 
0. In the limit of large Hilbert space dimensions (small h), for parameters such that 
the classical description of the dynamics shows global chaos, the eigenstates and eigen¬ 
values of the quantum dynamics have the statistical properties of an ensemble of random 
matrices. The appropriate ensemble depends on the properties of the quantum system 
under time-reversal 0. The ensemble of random matrices used to describe the Hamil¬ 
tonians unrestricted by the time reversal symmetry is the Gaussian Unitary Ensemble 
(GUE). Similarly, the ensemble of random matrices used to describe the Hamiltonians 
having a time reversal symmetry are given by the Gaussian Orthogonal Ensemble (GOE). 
The other class of random matrices typically studied are the random unitary matrices. 
They are employed for periodically driven systems, as models of the unitary “Eloquet” 
operators, F, describing the change of the quantum state during one cycle of the driving. 
Powers of the “Eloquet” operator, F”, give us a stroboscopic description of the dynamics. 
The ensemble of random unitaries are also known as the “circular ensembles”, originally 
introduced by Dyson 0. As was the case for random Hermitian matrices, time rever¬ 
sal symmetry arguments play a similar role in the choice of the appropriate ensemble of 
random unitaries employed to model the “Eloquet” operator to study the properties of the 
chaotic system. Depending on whether the system has time reversal symmetry or not, 
the appropriate ensemble of random unitaries is called the Circular Orthogonal Ensem¬ 
ble (COE) or the Circular Unitary Ensemble (CUE) respectively. The eigenvectors of the 
COE and CUE have the same properties as that for the respective GOE and GUE, but the 
eigenvalues are distributed differently. The circular unitary ensemble (CUE) is just the 
ensemble of random unitary matrices picked from U (n) according to the Haar measure. 
CUE eigenvalues lie on the unit circle in the complex plane, and hence the name. 
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From this fact, an important quantum signature of chaos was obtained by Bohigas and 
collaborators 0, describing the spectral statistics of quantum Hamiltonians whose clas¬ 
sical counterparts exhibit complete chaos using random matrix theory. Such signatures 
of quantum chaos have mainly focussed on the time-independent Schrodinger’s equation 
and features like energy spectra and eigenstates. 

Though quantum systems show no exponential separation under the evolution of a 
known unitary evolution, they do show a sensitivity to the parameters in the Hamilto¬ 
nian a. Peres Q showed that the evolution of a quantum state is altered when a small 
perturbation is added to the Hamiltonian. As time progresses, the overlap of the per¬ 
turbed and unperturbed states gives an indication of the stability of quantum motion. It 
was shown that if a quantum system has a classically chaotic analog, this overlap has a 
very small value. On the other hand, if the classical analog is regular, the overlap remains 
appreciable. In another perspective, as seen in the work of Schack and Caves, quantum 
systems exhibit chaos when they are perturbed by the environment. They become hy¬ 
persensitive to perturbations i), as seen in the information-theoretic studies of the cost to 
maintain low entropy in the face of loss of information to the environment. This particular 
feature of quantum chaotic systems has several interesting consequences. For example, 
Shepelyansky has done extensive work on the issue of many-body quantum chaos in the 
quantum computer hardware and its effect on the accuracy of quantum computation a 
in the absence of error correction. Recently, classical simulations of quantum dynamics 
have been connected to integrability and chaos Go). 

It is imperative to mention the role played by quantum information theory in the above 
journey. Quantum information science has added a whole new perspective to the study 
of quantum mechanics. This has resulted in a better understanding of quantum phenom¬ 
ena like entanglement and decoherence, and given us the tools to view certain quantum 
properties of physical systems as a resource. This has also enabled us to address the 
key questions in quantum chaos from a new perspective. As mentioned above, this has 
led to an information theoretic characterization of quantum chaos [|8] and explained the 
exploration of the behavior of chaotic quantum systems in the presence of environment 
induced decoherence in along with its connection to the quantum-to-classical transi¬ 
tion. The study of quantum chaos from a quantum information perspective is also closely 
related to the theory and application of random quantum circuits m. In the last two 
decades, quantum information theory has given us a new perspective in finding the fin¬ 
gerprints of chaos in quantum mechanics. The dynamical generation of entanglement and 
discord and information gain in tomography have been studied as signatures of classical 
chaos in the quantum world EHia. 

In ll24l . it was shown that information gain about an initial quantum state in the process 
of quantum tomography is a metric to characterize and quantify quantum chaos. In this 
work, we review this new information-theoretic characterization of chaos and show how 
this procedure can be used to distinguish between symmetry classes of various quantum 
maps. 

Quantum tomography is the process of estimating an unknown quantum state from the 
statistics of measurements made on many copies of the state. In this work, we extend 
our efforts on information gain in quantum tomography to characterize the properties of 
the underlying dynamics. In particular, we give new analytical results for the information 
gain for different classes of quantum maps depending on their time reversal and par¬ 
ity symmetry properties. The standard way to perform quantum tomography is to make 
projective measurements of an “informationally complete” set of observables and repeat 
them many times. The statistics obtained are used to estimate the expectation values of 
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the observables and hence the unknown initial state. 

Projective measurements pose a hurdle in exploring the connections between informa¬ 
tion gain in tomography and chaos due to large measurement back-action on the system. 
However, we overcome this by employing the protocol for tomography via weak continu¬ 
ous measurement developed by Silberfarb et al. ||25]. In this protocol, the ensemble is col¬ 
lectively controlled and probed in a time dependent manner to obtain an “informationally 
complete” continuous measurement record. We consider the case of a very weak mea¬ 
surement such that the back-action is negligible. This is possible when the uncertainty 
in any measurement outcome is small compared to the quantum uncertainty associated 
with the probe itself. We accurately model all of the quantum dynamics occurring in the 
system, and then use the measurement time history to give us information about the initial 
quantum state. The dynamics is “informationally complete” if the time history contains 
information about an arbitrary initial condition. Our goal is to characterize and quantify 
the performance of tomography, when the dynamics driving the system are chaotic in the 
classical limit. We use this to draw connections between the role played by regular and 
chaotic dynamics as well as the nature of symmetries of the dynamics in the tomography 
procedure. The work presented in this paper is intimately related to the protocols that 
have recently been implemented in the laboratory ESI . 

The remainder of this paper is organized as follows. In Sec. 2, we review the protocol 
for tomography via weak continuous measurement developed by Silberfarb et al. Il25l and 
Riofrio et al. Il27l . In Section 3, we demonstrate how information gain while performing 
tomography is a quantum signature of classical chaos. We perform numerical simulations 
of the reconstruction hdelity and its relationship to the degree of chaos in the dynamics 
that drive the system. We also show how the hdelity obtained and the corresponding 
metrics to quantify information gain can be used to distinguish quantum maps belonging 
to different symmetry classes. We then explain these results in terms of the properties of 
random states in Hilbert space. Our results are discussed and summarized in Sec. 4. 


2. Tomography via Weak Continuous Measurement 

In this section, we review tomography via a continuous measurement protocol. Consider 
an ensemble of N, noninteracting, simultaneously prepared quantum systems in an iden¬ 
tical, but unknown, state described by the density matrix pQ. Our goal is to determine po 
by continuously measuring an observable Oq. The ensemble is collectively controlled and 
probed in a time-dependent manner to obtain an “informationally complete” continuous 
measurement record. In order to achieve informational completeness, when viewed in the 
Heisenberg picture, the set of measured observables should span an operator basis for po- 
For a Hilbert space of hnite dimension d, and hxing the normalization of po, the set of 
Hermitian operators must form a basis of su{d). The measurement record is inverted to 
get an estimate of the unknown state. Laboratory realization of such a record is intimately 
tied to controllability, i.e., designing the system evolution is such a way as to generate 
arbitrary unitary maps. While it is desirable to obtain an informationally complete mea¬ 
surement record, we shall see that we can obtain high fidelity in tomography in some 
cases even when this is not the case ll28l . 

In an idealized form, the probe performs a QND measurement that couples uniformly to 
the collective variable across the ensemble and measures J2n=j ^ strong QND 

measurement, quantum backaction will result in substantial entanglement among the par¬ 
ticles. For a sufficiently weak measurement, the noise on the detector (e.g., shot noise of 
a laser probe) dominates the quantum fluctuation intrinsic to the measurement outcomes 


4 


of the state (projection noise). In this case, we can neglect the backaction on the quantum 
state and the ensemble remains factorized. In order to obtain a measurement record that 
can be inverted to reconstruct an estimate of the initial state, one must drive the system 
by a carefully designed dynamical evolution that continually maps new information onto 
the measured observable. In order to do so, the system is manipulated by external fields. 
The Hamiltonian of the system, H{t) = H[(j)i{t)], is a functional of a set of time depen¬ 
dent control functions, so that the dynamics produces an informationally complete 

measurement record A4. 

Then we can write the measurement record obtained as 


M{t)=Tr{Oop{t))+aW{t), (1) 

amplified by the total number of copies (N atoms in this case). Here aW(t) is a Gaussian- 
random variable with zero mean and variance cr^, which accounts for the noise on the 
detector. Since our goal is to estimate the initial state from the measurement record and 
the system dynamics, we will work in the Heisenberg picture. Rewriting Eq. [^in the 
Heisenberg picture, we get 

M(t)=Tr(0(t)po) + aW(t). (2) 

We sample the measurement record at discreet times so that 

M^ = Tr(0,po)+^W^. (3) 

Thus, the problem of state estimation is reduced to a linear stochastic estimation problem. 

The goal is to determine po, given {Mi] for a well chosen {Oi}, in the presence of 
noise {Wi}. We use a simple linear parametrization of the density matrix 


j tl — J. 

Po = ^ I'aEa, (4) 

Q —1 

where d is the dimension of the Hilbert space, are d^ — 1 real numbers (the components 
of a generalized Bloch vector), and [Ea] is an orthonormal Hermitian basis of traceless 
operators. We can then write Eq. [^as 

d^-i 

raTr(0,E^) + aW,, (5) 

a— 1 

or, in the matrix form as 


M = dr + crW, (6) 

which in general is an overdetermined set of linear equations with d^ — 1 unknowns 
r = (ri, ...,rd 2 _i). 

The conditional probablity distribution for the random variable M, given the state r, is 
the Gaussian distribution 

75(M|r) cx exp(-^(M - dr)^(M - dr)). (7) 

2 ( 7 ^ 
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We can use the fact that the argument of the exponent in Eq. |^is a quadratic function of 
r to write the likelihood function (ignoring any priors) as 

V{r\M) oc exp{-^{r - VmlV{ r - Tml)), (8) 

a Gaussian function over the possible states r centered around the most likely state, Vml, 
with the covariance matrix given by C = The uncontrained maximum 

liklihood solution is given by 

Vml = {d^dy^d^M. (9) 

The measurement record is informationally complete when the covariance matrix has full 
rank, — 1. If the measurement record is incomplete and the covariance matrix is not 
full rank, we replace the inverse in Eq. |^with the Moore-Penrose pseudo inverse ll29l . 
The eigenvectors of represent the orthogonal directions in operator space that we 
have measured up to the hnal time, and the eigenvalues determine the uncertainty, or the 
signal-to-noise ratio, associated with those measurement directions. 

When we have an incomplete measurement record, or in the presence of noise, the 
unconstrained maximum likelihood procedure does not give a density matrix that corre¬ 
sponds to a physical state. The estimated density matrix might have negative eigenvalues. 
We correct this by hnding a valid density matrix that is “closest” to pml, the density 
matrix obtained by the unconstrained maximum likelihood procedure. 

3. Information Gain in Tomography 

3.1 Metrics to quantify information gain 

Our protocol for quantum tomography via continuous measurement of a driven system 
gives us a window into the complexity of quantum dynamics and its relationship to 
chaos. Moreover, the experimental implementation of tomography by continuous mea¬ 
surement provides a useful platform for exploring these ideas in the laboratory ll26l . 
Quantum tomography deals with the extraction of information about an unknown quan¬ 
tum state through measurements. In our attempt to study chaos under this paradigm, we 
dehne metrics to quantify this information gain. These metrics characterize the ability 
of our control dynamics to generate a sufficiently high signal-to-noise ratio for measure¬ 
ments in different directions of the operator space. As we shall see, these metrics elucidate 
the connection between the degree of chaos and the fidelities obtained in tomography. We 
can quantify the information gain in a number of ways. 

1) Fidelity of Tomography: Eidelity of the reconstruction obtained in tomography is 
a metric for information gain which determines the degree of closeness of quantum states 
and is intimately related to how much information is obtained during the process. The 
fidelity is simply given by the overlap of the initial and the reconstructed state vectors. 
The fidelity between a target pure state |^) and the reconstructed state pis F = {y p \ ij}). 

2) Fisher Information (FI) of the Measurement Record: We can further quantify 
the correlation between chaos and the performance of quantum state estimation using 
information-theoretic metrics. The information obtained in measurement of a quantum 
system can be expressed in terms of the uncertainty of the outcomes summed over a set 
of mutually complementary experiments 130| . In terms of the Hilbert-Schmidt distance 
between the true and estimated state in quantum state reconstruction, averaged over many 
runs of the estimator, this information can be written as / = {Tr{{po — p)^}) OTI . which 
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in terms of the total uncertainty in the Bloch vector components is / = 

The Cramer-Rao bound tells us that this uncertainty obeys 

((Ar„)2)> (10) 

where F is the Fisher information matrix associated with the conditional probability dis¬ 
tribution, Eq. 1^ and thus I > Tr 

In general, for a multivariate parameter estimation problem, the Cramer-Rao bound 
gives 

Cove(T(A)) >F-\ (11) 

where the matrix inequality, A > B, 'k understood to mean that the matrix, A — B, 'k 
positive semidefinite. Here A is a d-dimensional random vector that contains information 
about the multivariate parameter, 9 = [0i,02, T(X) is the unbiased estimator of 

the multivariate parameter, and Cove(T(A)) is the covariance matrix of a set of unbiased 
estimators for the parameters 9. It quantifies the error in our estimation process. F is the 
multivariate generalization of the FI lf32l . 

Fmn = E(^—log/(a:;6»)—log/(a;;6»)), (12) 

C/Uyn 

where f(x] 9) is the probability density of the random variable X conditioned on the value 
of 9, and E denotes the expectation value. In the limit of negligible quantum backaction, 
we saturate this bound. This is because our probability distribution is Gaussian, regardless 
of the state. In that case, the Fisher information matrix equals the inverse of the covariance 
matrix, F = C“^, in units of and thus the Cramer-Rao bound reads 

Cove(T(A)) > C. (13) 

We consider the basis in which F, and hence C~^, is diagonal, 

f' = UFU^. (14) 

Such a transformation is provided by U composed from the eigenvectors of C. In this 
representation, the estimate of the newly transformed parameters fluctuate independently 
of each other. This suggests the possibility to form a single number that quantifies the 
performance of the tomography scheme as a whole by adding those independent errors, 
e, as 

e>Tr{C). (15) 

Thus, which is the collective FI, serves as a measure of the amount of information 

about the parameter 9 that is present in the data. 

3) Shannon Entropy of Eigenvalues of the Inverse Covariance Matrix: The mutual 
information, I[r; M], quantifies the information we have about parameters r from mea¬ 
surement record M, which is given by I[r; M] = iT(M) — FI (M|r) ||3^ . Here iJ is the 
Shannon entropy of the given probability distribution. The entropy of the measurement 
record, iF(M), arises solely due to the shot noise in the probe, and hence is a constant. 
The mutual information between the Bloch vector and a given measurement record can 
be expressed as the entropy of the conditional probability distribution (Eq. |7]l 

I[r;M] = -H{M\r) = -ilog(detC) = log(l/F), (16) 

where V is the volume of the error-ellipsoid whose semi-major axes are defined by the 
covariance matrix. 
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3.2 The Quantum Kicked Top 

How does the presence of chaos in the control dynamics influence our ability to perform 
tomography? In order to address this question, we chose the “kicked top” dynamics as 
the paradigm to explore quantum chaos in tomography. The Hamiltonian for the kicked 
top (after setting Ti= 1) is given by 


1 1 ^ 

H{t) = -pJa: + —kJ^ S{t-nT). 
r 2? 

n— — oo 


(17) 


Here, the operators, J^, Jy and are the angular momentum operators obeying the 
commutation relation [Ji, Jj] = iCijkJk The first term in the Hamiltonian describes a 
precession around the x axis with an angular frequency and the second term describes a 
periodic sequence of kicks separated by time period r. Each kick is an impulsive rotation 
about the z axis by an amount proportional to J^- Choosing the external field to act in 
delta kicks allows us to express the Floquet map (transformation after one period) in a 
simple form of sequential rotations as 


Ur = e 


ZOtJX 


( 18 ) 


where a and A are related to p and k, respectively, in terms of the kicking period. The 
evolution of the initial quantum state has the form where n enumerates the kick 

number or the periodic application of the map. The classical map can be obtained by 
considering the Heisenberg evolution of the expectation values of the angular momentum 
operators in a familiar way 0 . 

The classical dynamics consists of the motion of a unit spin vector on the surface of 
the sphere. The z-component of a spin and the angle (j), denoting its orientation in the x-y 
plane, are canonically conjugate, and thus the spin constitutes one canonical degree of 
freedom. The classical dynamical map has the same physical action as described above in 
the quantum context - precession of the spin around the x axis with an angular frequency 
a followed by an impulsive rotation around the z axis by an amount proportional to 
with a proportionality constant A. In our analysis, we fix a = 1.4 and choose A to be our 
chaoticity parameter. As we vary A from 0 to 7, the dynamics change from highly regular 
to completely chaotic. Since the total magnitude of the spin is a constant of motion, 
our classical map is two dimensional. We visualize the phase by plotting the z and y 
components of motion after every application of the dynamical map. 

Figure [T] shows four different regimes of classical dynamics. With the parameters a = 
1.4, A = 0.5 (Fig. [T^), the dynamics are highly regular. When a = 1.4, and A = 2.5 
(Fig. [^), we see a mixed space with chaotic and regular regions of comparable size. The 
parameters, a = 1.4, A = 3.0 (Fig. [^), give a phase space that has mostly chaotic regions 
and finally, a = 1.4, A = 7.0 gives a completely chaotic phase space (Fig. [^). 

A central result of quantum chaos is the connection with the theory of random matrices 
0 . In the limit of large Hilbert space dimensions (small h), for parameters such that the 
classical description of the dynamics shows global chaos, the eigenstates and eigenvalues 
of the quantum dynamics have the statistical properties of an ensemble of random matrices 
0 . The appropriate ensemble depends on the properties of the quantum system under 
time-reversal symmetry 0. We thus seek to determine whether there exists an anti-unitary 
(time reversal) operator T that has the following action on the Floquet operator, 

TUrT-^ = Ul= ( 19 ) 
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Figure 1. Phase space plots for the kicked top in four regimes, (a) Regular phase 
space: a = 1.4, A = 0.5, (b) Mixed phase space: a = 1.4, A = 2.5, (c) Mostly 
chaotic: a = 1.4, A = 3.0, (d) Fully chaotic phase space: a = 1.4, A = 7.0. The 
figures depict trajectories on the southern hemisphere (x < 0) of the unit sphere where 
X = —,Y — ^ and Z = —, and we take the limit j —>■ oo to get the classical limit 
as in |4|. 
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Considering the generalized time reversal operation 


T = 


( 20 ) 


where K is the complex conjugation operator. It then follows that 


TUrT-'^ = 





( 21 ) 



e 


—iaJx 


so the dynamics is time-reversal invariant. Moreover as = 1, there is no Kramer’s 


degeneracy. Given these facts, for parameters in which the classical dynamics is globally 
chaotic, we expect the Floquet operator to have the statistical properties of a random 
matrix chosen from the circular orthogonal ensemble (COE) El- 

In order to have maximum information gain, we need to condition the dynamics so that 
we maximize 1/V = -y/det (C“^). The quantity Tr{C~^) is constrained at f„. One can 
show that after n steps 



( 22 ) 


where ||(!I(0)|p = is the Hilbert-Schmidt square norm, and 0(0) = 

Jz for our case. Therefore, from the theorem of the arithmetic and geometric means. 



(23) 


where U = — 1 is the rank of the regularized covariance matrix. The maximum 

possible value of the mutual information is attained when all eigenvalues are equal, sat¬ 
urating the above inequality. At a given time step, the dynamics that gives the largest 
mutual information is the one that makes the eigenvalues most equal. If we normalize the 
eigenvalues of the inverse of the covariance matrix, then as a probability distribution, its 
Shannon entropy E, is a measure of how evenly we have sampled all the directions in the 
operator space. We reach the maximum entropy when we have measured all directions in 
the space of matrices equally, Emax = log((i^ — 1). This is the most unbiased measure¬ 
ment we can implement that will lead to the highest fidelities, on average, for a random 
state. 

The collective FI, rj. , tells us about the amount of information our measurement 
record contains about the parameters that define the density matrix. Figure ffl shows the 
behavior of the FI as a function of the number of applications of the kicked top map, and 
for different values of the chaoticity parameter. We see that the rate of increase of the FI 
is correlated with the degree of chaos present in the control dynamics. As our dynamics 
become increasingly chaotic, we obtain higher values for the FI at a given time. We expect 
the FI to be correlated with the average fidelities of estimation for an ensemble of random 
states. 
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Fidelity 



Figure 2. Fidelity of reconstruction as a function of the number of applications of the 
kicked top map. The fidelity is calculated as the average fidelity of reconstruction of 
100 states picked at random according to the Haar measure. The parameters of the 
kicked top are as described in the text, with a = 1.4 fixed. We show the fidelity for 
different choices of the chaoticity parameter. Both the rate of growth and the final value 
of the fidelities are increased with higher values of A. 
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Figure 3. The Shannon entropy of the normalized eigenvalues of the inverse of covari¬ 
ance matrix as a function of the number of applications of the kicked top map: (a) Short 
time behavior (b) long time/asymptotic behavior. The parameters are as described in 
the text. 



Figure 4. The FI of the parameter estimation in tomography as a function of the 
number of applications of the kicked top map: (a) Short time behavior (b) long 
time/asymptotic behavior. The parameters are as described in the text. 
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3.3 Signatures of chaos : Information gain in the fully chaotic regime and random ma¬ 
trix theory 

3.3.1 Results and Discussion 

We are now ready to explore the role of chaos in the performance of tomography. Through¬ 
out this section, we consider spin J = 10, a d = 21 dimensional Hilbert space, which is 
sufficiently large that a minimum uncertainty spin coherent state is a sufficiently confined 
“wavepacket” that it can resolve features in the classical phase space. Figure|^shows the 
average hdelity of reconstruction of 100 states picked at random according to the Haar 
measure as a function of the number of applications of the kicked top map, and for differ¬ 
ent values of the chaoticity parameter. We see that the rate of increase in hdelity increases 
with the degree of chaos. The hnal hdelity achieved after a hxed number of kicks is also 
correlated with the degree of chaos. We can understand the above results by studying the 
information gain in tomography as a function of the degree of chaos in the control dynam¬ 
ics. Figure 1^ shows the behavior of the entropy E of the covariance matrix, as dehned 
above, as a function of the number of applications of the kicked top map, and for different 
values of the chaoticity parameter. We see that the rate of increase of entropy for short 
times. Fig. Si, is correlated with the degree of chaos present in the control dynamics. 
The asymptotic value of the entropy reached also increases with the chaoticity parameter. 
Chaotic dynamics provides a measurement record with a large signal-to-noise ratio in all 
the directions in the operator space. An increase in the chaoticity parameter results in an 
increasingly unbiased measurement process that will yield high hdelities for estimating 
random quantum states. Figure]^ shows the behavior of the entropy at short time scales, 
while we see asymptotic behavior in Fig. [^. 

When the system is driven by dynamics that are completely chaotic, we expect the 
information gain and the hdelity to follow the predictions from random matrix theory. 
Figure]^ shows the behavior of the hdelity. Shannon entropy and the FI of the inverse of 
the covariance matrix as a function of the number of applications of the kicked top map 
(the blue line) and compares them with the corresponding quantities for a typical random 
unitary picked from the COE (the green line). We see a strong agreement between our 
predictions from random matrix theory and the entropy calculation for the evolution by a 
completely chaotic map. 

We test our predictions from the random matrix theory for chaotic maps without a 
time reversal symmetry. For example, another type of the “kicked top” map without time 
reversal symmetry lf33l is given by 

JJ — — — (24) 

In Fig. 1^ we repeat the above calculations for this map. In this case, the appropriate 
random matrix ensemble is the CUE. We see an excellent agreement between the behavior 
of the hdelity. Shannon entropy and the El, as predicted by random matrix theory, and that 
for the evolution by a completely chaotic map without the time reversal symmetry ll^ . 

When all the eigenvalues of the inverse of the covariance matrix are equal, we have an 
upper bound on the entropy, Emax = log((i^ — 1). Eigures|^-|^ compare the entropy 
values achieved by the repeated application of the same unitary (time reversal invariant or 
otherwise) to Emax- We see that we fall signihcantly short of Emax by such a procedure. 

So far, we have considered the application of the same unitary matrix periodically to 
obtain the measurement record. However, this alone does not give us an informationally 
complete measurement record; high hdelities are reached only when we make use of 
the positivity constraint. On the other hand, we can consider application of a series of 
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t„ (a) t„ (b) t. (c) 

Figure 5. Comparison between the tomography performed by the repeated application 
of kicked top in the fully chaotic regime (the blue line) and that by a typical random 
unitary picked from the COE (the green line), (a) The average fidelity of reconstruc¬ 
tion of 100 states picked at random according to the Haar measure, (b) The Shannon 
entropy of the normalized eigenvalues of the inverse of covariance matrix as a function 
of the number of applications of the map. The dotted line gives the upper bound on the 
entropy, Emax = log(rf^ - 1). 




Figure 6. Same as Figure 4, but for a kicked top without time reversal invariance (Eq. 
\24\ (the blue line). In this case as well, the results are well predicted by modeling the 
dynamics by random matrices sampled from the CUE (the green line). 
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different random unitaries 0^ . In that case, we expect to rapidly reach an informationally 
complete set and thus rapidly gain information about tomography. In Fig. we plot the 
fidelities, Shannon entropy and the FI achieved by applying a different random unitary, 
picked from the unitarily invariant Haar measure, and compare it with the results obtained 
by the repeated application of the same unitary (picked from the COE and CUE). We also 
see that we reach the upper bound, Emax, asymptotically, by this method. Indeed, an 
application of a different random unitary is the most unbiased dynamics we can hope to 
perform. 





Figure 7. Comparison between tomography performed by applying a different random 
unitary at each time step, picked from the unitarily invariant Haar measure (magenta 
line) and that by a repeated application of a random unitary picked from the COE (blue 
line) and the CUE (red line) . The dotted line gives the upper bound on the entropy, 
Emax — log(d' - 1) 


3.3.2 Analytical expressions for Information gain 

In this section, we use random matrix theory to predict the information gain in tomogra¬ 
phy when we apply the unitary map Ur periodically. 

1) The quantum kicked top : In our system, we have a symmetry given by the parity 
operator R that has the form 

R = . (25) 

Our unitary map, U, being the kicked top or the appropriate COE sampled matrix, will 
commute with R, i.e., [R, U] = 0. Thus, there exists a basis in which U and R are 
diagonal. Eirst, note that the eigenvalues of R are ±1. Then, let’s define a basis, {|i?y)}, 
where \Rj) = \Rj~^} for j = l,...,a, and \Rj) = \Rj'^^) for fc = a + 1,. ■ ■ ,d, 
corresponding to the eigenvalues —1 and -fl, and where a G {(d + l)/2, (d — l)/2}. 
Since U is also diagonal in this basis, an asymptotic approximation to the inverse of the 
covariance matrix is 


« n 


^ l(RklOolRfl^lRk,Rj)(Rk,E, 


j,k-l 


+ (Ejl Oo jRj) (Rkl Oo jRk) iRj, Rj) (Rk, E/cl 

j^k=l 


( 26 ) 
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Our initial observable, Oo = Jz, anti-commutes with R, meaning that 

= _ j^. 


(27) 


Because of this, we see that 

= 0, forj,k = l,...,a, (28) 

and 

{R^+^\Jz\Ri^^) = -{Rl+^\Jz\Ri^'>) =0, forj,fc = a + l,...,d. (29) 

As discussed above in Eq. we know that after time f„, the trace of the inverse covari¬ 
ance matrix is given by Tr(C“^) = n/3, where /3 = |jO(0)|p is a constant independent 
of the Floquet map, Ur, driving the system. 

Thus, the matrix representation of Jz in the ordered basis in which R is diagonal is 
anti-block diagonal. We immediately see that Eq. [26]simplihes to 




Y,mk\Oo\R,)?\Rk,Rj) {Rk,R, 


(30) 


In this basis, C~^ is approximately diagonal. So we can actually give an analytical for¬ 
mula for the Shannon entropy. Remember that we previously dehned the normalization 
factor as fi. So the eigenvalues of C~^ are simply |(i?fe| Oq |7?j)P//3. To compute the 
expected value of the Shannon entropy, we use the results of Wootters llTSll for the ex¬ 
pected value of entropy of the entries of a state expressed in a random basis, and sampled 
from the appropriate ensemble. We see that since Jz is anti block diagonal, there are only 
2 X (d — l)/2 X (d -f l)/2 nonzero terms. 

Now, we can directly use Wootters formula for the expected Shannon entropy, 

Hexp = log(i9) - 0.729637 = log ^ - 0.729637. (31) 


For d = 21, we get idexp = 4.66. Numerically, we hnd a somewhat larger value for 
the kicked top, Hkt = 4.85 and Hav = 4.69 for entropy averaged over 100 block 
diagonal COE matrices. This is due to the fluctuations in H about the expected value and 
these fluctuations reduce as we increase d and we And an excellent convergence with the 
analytical expression derived above. 

2) The CUE ; The above analysis can be carried over to a quantum map that does not 
have time reversal symmetry. In this case, since there is no parity symmetry, there are 
d^ — 1 nonzero terms in Eq. 26 and therefore, we get for the expected Shannon entropy 
as 


Hexp = log {d? - 1) - 0.729637, (32) 

which agrees remarkably well with our numerical simulations (for d = 21, Hexp = 5.35). 

3) A different Haar random unitary at each time step; In this case, we explore the 
complete Hilbert Space and we get Hexp = log(d^ — 1), which agree very well with 
our simulations (for d = 21, Hexp = 6.08). The maximum possible value of the mutual 
information is attained when all eigenvalues of the covariance matrix are equal. In order to 
extract the maximum information about a random state, we must measure all components 
of the Bloch vector with maximum precision. In finite time, we obtain the best estimate 
by dividing evenly among all observables. 
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4. Conclusion and Outlook 


The missing information in deterministic chaos is the initial condition. A time history 
of a trajectory at discrete times is an archive of information about the initial conditions 
given perfect knowledge about the dynamics. Moreover, if the dynamics is chaotic the 
rate at which we learn information increases due to the rapid Lyapunov divergence of 
distinguishable trajectories and we expect unbiased information because of the ergodic 
mixing of phase space. That is, if the information is generated by chaotic dynamics, the 
trajectory is random, and all initial conditions are equally likely until we invert the data 
and discover the initial state. 

Dynamics sensitive to the initial conditions will reveal more information about the ini¬ 
tial conditions as one observes the system trajectory in the course of time. Classically 
chaotic dynamics generates this unpredictability, or information to be gained about the 
initial coordinates of the trajectory. Similarly, we found that the rate at which one obtains 
information about an initially unknown quantum state in quantum tomography is corre¬ 
lated with the extent of chaos in the system. This is a new quantum signature of classical 
chaos. In fact, our results can be regarded as signatures of chaos in quantum systems un¬ 
dergoing unitary evolution, as measurement backaction is negligible. We have been able 
to quantify the information gain using the FI associated with estimating the parameters 
of the unknown quantum state. When the system is fully chaotic, the rate of information 
gain agrees with the predictions of random matrix theory. 

At its core, our approach is akin to the Kolmogorov-Sinai (KS) entropy measure of 
chaos ll^ . Incomplete information about the initial condition leads to unpredictability of 
a time history. In the presence of classical chaos, in order to predict which coarse-grained 
cell in phase space a trajectory will land at a later time, we require an exponentially 
increasing fine-grained knowledge of the initial condition. The KS entropy is the rate of 
increase, and is related to the positive Lyapunov exponents of the system. Is there a mean¬ 
ingful quantum definition of KS entropy? Our results seem to suggest this. In order to 
predict the measurement record with a fixed uncertainty, we need to learn more and more 
about the initial condition. Is the rate at which we obtain this information exponentially 
fast when the system is quantum chaotic? Does this converge to the classical Lyapunov 
exponents in the limit of large action (small Ti)l There are many important subtleties in 
these questions. 

As we gain more and more information, eventually quantum backaction becomes im¬ 
portant in the measurement history. The number of copies we have and the shot noise 
on the probe limits the ultimate resolution with which we can deduce the quantum state 
Ell. Unlike classical dynamics, we can never consider infinite resolution, even in princi¬ 
ple. The quantum resolution is limited by the size of Ti. As the dimension of the Hilbert 
space increases, and hence the effective h decreases, we expect to see an even sharper 
difference in the information gain as a function of chaoticity. In the limit when d, the 
dimension of the Hilbert space, becomes infinity, we expect the rate of information gain 
to be intimately related to the classical Lyapunov exponents. How all this translates into 
a quantum definition of KS entropy is an important subject of further investigation. 

In principle we never have perfect knowledge of the dynamics. This is related to hyper¬ 
sensitivity to perturbations ll38l in quantum chaotic dynamics. This implies that, though 
quantum systems show no sensitivity to initial conditions, due to unitarity, they do show 
a sensitivity to parameters in the Hamiltonian llsl l^ . How does this fundamentally limit 
our ability to perform quantum state reconstruction when the system is sufficiently com¬ 
plex, and the equivalent dynamics is chaotic. This poses interesting questions for quantum 
tomography and, more interestingly, for quantum simulations. Under what conditions are 
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the system dynamics sensitive to perturbations and how does this effect our ability to per¬ 
form quantum tomography? Under what conditions does the underlying quantum chaos 
affect our ability to accomplish quantum simulations in general? We hope to address these 
questions in our future work. 
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